Index of refraction of gases for matter waves: 
effect of the motion of the gas particles on the calculation of the index. 
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Two different formulae relating the index of refraction n of gases for atom waves to the scattering 
amplitude have been published. We show here that these two formulae are not consistent with the 
definition of the total scattering cross-section while the formula developed by one of us (C.C.) in 
her thesis is in agreement with this standard knowledge. We discuss this result, in particular in 
the neutron case for which such an index was first introduced. We finally evaluate the index of 
refraction as a function of well known quantities and we discuss the order of magnitude of the ratio 
of (n — l)/nt, where nt is the gas density. 



I. INTRODUCTION 

To describe the propagation of matter waves through 
a medium, it is possible to use an index of refraction, as 
commonly done for light. This idea was first introduced 
around 1940 for the propagation of neutrons through 
matter (see the review papers by L. L. Foldy [l[ and by 
M. Lax 0). 

Atom intcrfcromctry has permitted to study the prop- 
agation of an atom wave through a dilute gas and the 
first measurements of the index of refraction of gases for 
sodium waves were made in 1993 at MIT by J. Schmied- 
mayer et al. [j| H| who measured the attenuation and the 
phase shift of the transmitted wave. Further works in the 
same laboratory have led to the observation of glory os- 
cillations of the index as a function of the sodium atom 
velocity @,[El3- More recently, our group has measured 
the index of refraction of several gases for lithium waves 

a. 

Whatever is the nature of the wave and of the medium, 
the index of refraction describes the modification of the 
propagation of an incident wave due to the waves scat- 
tered in the forward direction by the particles of the 
medium: the scattered waves interfere with the incident 
wave and modify its phase and its amplitude. The mod- 
ification of the phase induces a modification of the wave 
velocity, described by the real part of the index of refrac- 
tion, while the modification of the amplitude is described 
by its imaginary part. 

In practice, the index of refraction is proportional to 
the complex forward scattering amplitude 0, licj . The 
imaginary part of this amplitude is related to the total 
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cross section, which is traditionally measured by beam 
scattering experiments, whereas its real part can be mea- 
sured only by atom interferometry. This amplitude ex- 
hibits resonances, for a collision energy comparable to 
the potential well depth, and glory oscillations, for larger 
energies. Glory oscillations are the quantum consequence 
of the existence of an undeflected classical trajectory due 
to the compensation of attractive and repulsive forces 




For atom waves, the forward scattering amplitude can 
be calculated if the interaction potential between a par- 
ticle of the wave and a particle of the target gas is 
known. Several papers [H QJ, QJ, M, M, El 111, M, HI 
have discussed the theory of the index of refraction 
for atom waves and their equations have been used 
to com par e calculated values of the index of refraction 
[H [3 [H [H, [H HI, EI H with experimental results. 
A detailed comparison is possible only if the thermal mo- 
tion of the target gas is not neglected and we show here 
that the formulae [HI, [l6| introduced to take into account 
this motion disagree with collision physics. We propose 
a formula developed by one of us (C. C.) in her thesis 
[13]: this formula agrees with standard results of colli- 
sion physics as well as with a recent calculation based on 
quantum Boltzmann equation due to K. Hornberger [24j . 

In the present paper, we first recall the previous re- 
sults concerning the index of refraction. Then, we explain 
why these formulae are not in agreement with the Beer- 
Lambert law and we extract from this discussion a new 
formula giving the imaginary part of the index of refrac- 
tion, which is generalized to the real part of the index 
of refraction. We then discuss the differences between 
the old and new formulae. We verify that our formula is 
in agreement with the results concerning the index of re- 
fraction for slow neutrons. Finally, we compare the order 
of magnitude of the index of refraction of gases for light 
and atomic waves. In the cases where the index of refrac- 
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tion for atom waves has been measured, the values of the 
ratio (n — l)/n t , n t being the target gas density, are very 
close to the value of this ratio for the index of refraction 
for ordinary gases for visible light. This similarity is in 
fact a coincidence, without any physical meaning. 



IV. INDEX FOR FIXED SCATTERING 
CENTERS 

When a plane wave of wavevector k p enters in a 
medium made of target particles t, its propagation is 
modified, with a different wavevector k PjTO given by: 



II. ORIGIN OF THE PROBLEM 

When one compares light waves and atom waves inter- 
acting with an atomic (or molecular) target gas, there is 
a large difference, which is precisely related to the motion 
of the target particles. 

In the case of light, the photon velocity, almost equal 
to c in dilute matter, is considerably larger than the tar- 
get particle velocity. At the same time, the photon mo- 
mentum is usually considerably smaller than the target 
particle momentum. The atom motion induces Doppler 
effect and the photon momentum induces atom recoil. 
Because the velocity of light is usually so much larger 
than the one of the target particles, these two effects 
have small consequences on the index of refraction, if we 
except the frequency range close to a sharp resonance 
line. This simple idea remains true, even if the target 
particles move rapidly in the laboratory, because one can 
always study the problem in their rest frame. 

In the case of atomic waves, the velocity of an atom of 
the wave is usually comparable to the velocity of a target 
particle and, at the same time, the scattering properties 
depend rapidly on the relative velocity. Moreover, the 
momentum of an atom of the wave and the one of a 
target particle are usually of comparable magnitudes. In 
practice, it is absolutely necessary to take into account 
the target particle motion to make a realistic calculation 
of the index of refraction. 



III. NOTATIONS 



In the laboratory frame, a plane wave of wavevector 
k p describes the propagation in vacuum of a projectile p 



of mass m p and velocity v p : 



Tlk p = 77/ 



p v p 



(1) 



The projectile can be any massive particle described by 
quantum mechanics, a neutron, an electron, an atom or 
a molecule. This wave propagates through a gas of den- 
sity nt made of target particles t of mass mt and velocity 
v t , the target particles being also described by quantum 
mechanics. The wavevector k r describes the relative mo- 
tion of the projectile p in the center of mass frame for a 
collision with a target particle t: k r = fj,v r /h, where v r = 
v p — V( is the relative velocity and \i = m p mt / (jn p + m t ) 
is the reduced mass. The present notations will be used 
to write all previous formulae in order to facilitate com- 
parison. 



nk„ 



(2) 



where n is the index of refraction. If the medium is de- 
scribed by a random distribution of fixed scattering cen- 
ters, the index of refraction n is given by: 



7i = 1 + 2imt 



/(M 

V2 
V 



(3) 



where f(k p ) is the forward scattering amplitude of the 
wave scattered by one scattering center. f(k p ) has the 
dimension of a length and, as usual, the index n is a 
dimensionless quantity. A general derivation of this for- 
mula is found in the paper of L. Foldy [l[ (see also the 
review of M. Lax Q and references therein). Equation 
(|3"|) is the first order term of an expansion in powers of 
the gas density n t and it is a good approximation if the 
following conditions are fulfilled: 

• the mean distance n t between nearest scatter- 
ing centers must be considerably larger than the 
projectile de Broglie wavelength A p = 2n/k p , i.e. 



n t 1/3 Ap < 1. 



• the mean distance n t between scattering cen- 
ters must be considerably larger than the range of 
the interaction potential. Unfortunately, as real- 
istic atom-atom interaction potentials V(r) vanish 
only for an infinite distance r, the range has not a 
clear definition for such potentials. 

• the density n t of scattering centers is low enough 
so that the mean field correction is negligible i.e. 
(n — 1) -C 1. This last condition, which involves a 
condition on k p and f(k p ), depends on the interac- 
tion potential and collision energy. 

Practically, for thermal atoms waves and a target gas 
near room temperature, with relative velocities of the or- 
der of 10 3 m/s and target densities up to n t ~ 10 19 m -3 
used in the experiments 0, 0, H, H, 0, @] , the mean in- 
teratomic distance n t 1 ' 3 is larger than 5 x 10~ 7 m, the 



< 10 



-10 



and 



index of refraction is of the order of \n — 1| 
these three conditions are well fulfilled. 

In the early days of atom interferometry, the target 
particles were treated as fixed scattering centers and the 
index was given by equation ([3]): this was the case of 
the first paper Q dealing with the index of refraction 
for atom waves, written by J. Schmiedmayer et al. in 
1993 and of the paper written by one of us (J.V.) in 1994 
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[121 ]. As soon as experimental values [J] of the index 
of refraction of gases for atomic waves became available 
in 1995, it appeared necessary to take into account the 
target gas thermal motion. The next section presents the 
equations used by different groups. 

V. FORMULAE GIVING THE INDEX OF A 
GAS AT THERMAL EQUILIBRIUM 

Three research groups have been involved in the calcu- 
lation of the index of refraction of gas for atomic waves 
and have worked on the proper way to take into account 
the motion of the targets. For convenience and clarity, 
the formulae used by each group are presented separately. 

A. Publications of D. Pritchard and co-workers 

After a first paper where targets were considered at 
rest 0, this research group described in 1995 a set of 
measurements of the index of refraction of gases for 
sodium waves [4]. To interpret their experiment, these 
authors use the following equation: 

n = i + 2 ™ ( ZiM (4) 

where the thermal average is not explicitly discussed. 
The same equation is also used in the review paper [5j 
written by this group in 1997. A thermal average is taken 
into account in the paper by T. D. Hammond et al. @ 
published in 1997, with an index of refraction given by: 



C. Publications of our research group in Toulouse 

In 1995, E. Audouard et al. published a calculation 
of the index of refraction of argon gas for sodium waves 
p^ |. This work was the first one by our group taking 
into account the effect of thermal averaging. We have 
made an error in the algebra which was corrected in [l4| . 
The derivation of the thermal average formula was given 
in a following paper by C. Champenois et al. [16]. Our 
calculation was based on the Fizeau effect: this effect is 
a frame dragging effect well known for light, which has 
also been studied in the case of neutron matter waves 
dESSzi- We had written: 

]%,m = k p + ((ncM - l)k r ) (6) 

with the center of mass index of refraction ucm given 
by: 

n CM = l + 27m t 4^ (7) 

As above, the brackets () stand for the average over the 
velocity distribution of the target particles. We will not 
recall the final result [l6j corresponding to a Boltzmann 
distribution of the target gas, because this result is not 
correct. After discussion with the group of A. Dalgarno 
in 1998, whose results are presented above, we were con- 
vinced that our formula was wrong, in particular because 
our description of the Fizeau effect is not correct. In her 
thesis 17], Caroline Champenois derived a new formula 
(see equations V.78 and V.81): 



1 + 27T 



/ f(K) 



(5) 



where the brackets () mean the average over the velocity 
distribution of the target particles. 



B. Publications of A. Dalgarno and co-workers 



n = 1 + 2im t - 



MO 



(8) 



As the derivation of this formula was very involved, we 
will not reproduce it here but, in the next section, wc 
use a simpler argument to convince the reader that this 
formula is the right one. 



Equation ([5]) appears for the first time in the paper of 
R. C. Forrey et al. [l5|, written in 1995 and published in 
1996. In their paper, R. C. Forrey et al. [15j calculate the 
distribution of the relative wavevector k r and, for sake of 
completeness, we reproduce this calculation in Appendix 
A. The group of A. Dalgarno has published a series of 
papers on the index of refraction for atom waves (20l . [23j 
or for electron waves [l|| and this work was continued by 
R. C. Forrey and co-workers (22J. 

In 2002, R. C. Forrey et al. published a paper pj] 
entitled 'On the statistical averaging procedure for the 
refractive index of matter waves'. In this paper, they do 
not give a complete derivation of their formula but they 
state that a key step is the Lorentz invariance of the ratio 
f(k r )/k r = f (kp) / kp. 



VI. DISAGREEMENT WITH CLASSIC 
RESULTS AND A NEW FORMULA 

In this part, we show that neither equation ([5]) pro- 
posed by R. C. Forrey et al. nor our equations © and 
([7|) are in agreement with well accepted atomic collision 
results. We first introduce the total scattering cross sec- 
tion and relate it to the imaginary part of the index of 
refraction. 



A. The Beer-Lambert law 

In nonrelativistic mechanics, the total scattering cross 
section a(v r ) is related to the number of collisions dN co u 
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occurring during a time dr in a volume dV between a 
projectile p and a target t: 



dN, 



coll 



drdV 



n p ntcr(vr)v r 



(9) 



From this equation, we deduce the mean number of col- 
lisions dN co u iP /dT encountered by a projectile p per unit 
time: 



dN, 



call, p. 



dr 



n t a(v r )v r 



(10) 



We now consider a beam of projectiles p crossing a slab 
of target gas, with the velocity v p perpendicular to the 
slab. A slab of thickness dL is crossed by a projectile p 
in a time dr = dL / v p and the mean number of collisions 
for a projectile is given by: 



dN, 



coll,p 



n t <j(v r ) — dL 



(ii) 



From this equation, one can deduce the transmission T 
of the slab i.e. the fraction of the incoming flux which 
has crossed the slab without any collision. For a finite 
thickness L, the transmission T is obtained by a straight- 
forward integration: 



with t = exp \—Tm(n — l)kL] and ip = lZe(n—l)kL. The 
transmission in intensity is T — t 2 , which depends solely 
on the imaginary part of (n— 1): 



T = exp [-21m(n - l)k p L] (16) 



C. Consequences for the index of refraction 

ASs the attenuation of the beam calculated in the two 
formalisms must be the same, the imaginary part of the 
index of refraction is related to the effective cross section 
by the following equation: 



Im(n — 1) = 



n t (a eff (v p )) 



(17) 



The total cross section is related by the optical theorem 
to forward scattering amplitude [101 ] : 



<r(v r ) = 4-7T 



Tm(/(Av)) 



(18) 



Using this relation and equations (fl~4|) and (jTTJ) , we ob- 
tain a formula giving the imaginary part of the index of 
refraction in agrement with Beer-Lambert law: 



T = exp 



-nt<j{v r )—L 



(12) 



This equation is the Beer-Lambert law, usually written 
with an effective cross section (T e ff(v p ) = a(v r )v r /v p : 



T = exp [-n t a e ff(v p )L] 



(13) 



If the target velocity is spread, with a normalized distri- 
bution P(v t ) (i.e. verifying J P(v t )d 3 v t — 1), we must 
replace the effective cross section a e ff(v p ) in equation 
(p~3|) by its average (& e ff( v p)) given by: 



(a eff (v p )} = I P(v t )aK)^d 3 v t . (14) 

Vn 



B. Wave-like description of the attenuation of a 
beam by a slab 



Xm(n — 1) = 27m t 



i P + mt I lm(f(k r j) ' 



m, 



V 



(19) 



where the brackets () mean the average over the target 
velocity distribution _P(v t ). Once we have an expression 
for the imaginary part of (n— 1), we get the real part by 
a simple generalization: 



1 + 27m t - 



mt //( /> , ) 



/■ 2 

p 



(20) 



This equation agrees with the result derived by C. Cham- 
penois in her thesis [ijj- It agrees also with a result re- 
cently obtained by K. Hornberger using the formalism 
of the quantum Boltzmann equation [24|. We want to 
point out that equation (|20[) can be applied with any 
type of velocity distributions, as could be produced by 
a gas flow or a molecular beam, and not only with a 
Maxwell-Boltzmann distribution corresponding to ther- 
mal equilibrium. 



We are going to calculate the transmission T of the 
beam through the same slab of length L, using the wave 
point of view. Let ipmc be the incident wave and iptrans 
the transmitted wave given by: 

^Ptrans — & ^ P ^inc — t& ^^Pinc (1^) 



D. Comparison of the new and previous formulae 

The new formula, equation (|20p . differs from equation 
© established by R. C. Forrey et al. [UGl: 

- the denominator is kf kip in equation §5§ and k 2 in 
equation (j2"D|) ; 
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- the mass ratio (m p + mt)/mt present in equation (1201) 
is absent from equation ([5]). 

- from equation ([5]) , we deduce an effective cross section 
given by: 



Weff(vp)) = / P{vt)a(v r )dv t 



(21) 



This expression of (c e //(u p )) differs from equation JTJ 
and in the v p <C a limit, where a is the thermal velocity 
defined in Appendix A, the effective cross section given 
by equation (|2"Tj) is independent of the projectile velocity 
V p , while the correct behavior given by equation (|14[) is 
(c e ff(vp)) oc 1/vp. This is a well known result, recog- 
nized in everyday's life: we run under the rain to get less 
wet! 

However, in the opposite limit when the target gas 
temperature vanishes, v< = so that v r = v p and k r = 
fiVp, it is easy to verify that equation (|20p and equation 
([5]) are then equivalent. This equivalence suggests that, 
even when the temperature does not vanish, the index 
calculated by these two formulae will not differ strongly 
as long as a C u p . 



VII. DISCUSSION OF THE NEUTRON CASE 



where Sq is the s-wave phase shift, a is the scattering 
length defined by a = — lim (tan<5o/£v) when k r — ► 0. If 
we keep only the leading term in equation (|22|) . f(k r ) « 
—a, our formula (|2H)) . is equivalent to the result of D. 
Klcinman and G. Snow 13 ill: 



1 — 2im t 



m„ + m t a 



rn L 



h 2 
V 



(23) 



As k r is absent from the result, there is no Doppler ef- 
fect on the index of refraction. However, if we take into 
account the first order term in k r in equation (|22p , the 
forward scattering amplitude has a non- vanishing imagi- 
nary part, Tm (f(k r )) « k r a 2 . As this imaginary part is 
linear in k r , the imaginary part of the index of refraction 
is sensitive to the motion of the target particles but this 
imaginary part, which is very small, is usually ignored. 
The fc r -dependence of the imaginary part of the index of 
refraction has no practical consequences but this remark 
proves that the absence of Doppler effect on the index of 
refraction of matter for neutrons is a very special prop- 
erty valid only for the real part of the index of refraction 
in the s-wave limit. 



VIII. COMPARISON OF THE ORDER OF 
MAGNITUDE OF THE INDEX OF REFRACTION 
OF GASES FOR LIGHT AND MATTER WAVES 



It is also interesting to apply equation (|2U)l to the well 
known case of neutron waves. Neutrons are scattered 
only by nuclei, if we except the case of ferromagnetic 
materials in which the magnetic interactions of the neu- 
tron spin cannot be neglected. In the low energy domain 
where the formalism of the index of refraction is useful, 
the neutron-nucleus scattering process is almost always 
dominated by s-wave scattering [H, [29[ . 

The index of refraction is frequently calculated as the 
consequence of an effective potential which is related to 
the scattering length of the neutron-nucleus interaction 
potential. However, there are few papers devoted to the 
theoretical relation between scattering theory and the in- 
dex of refraction for neutrons, in which the motion of the 
nuclei is taken into account. This is the case of the papers 
by B. A. Lippmann and J. Schwinger [30j and also of a 
brief note by D. Kleinman and G. Snow [31(, who state: 
"This derivation of the index clearly shows that there is 
no Doppler effect due to the motion of the nuclei, because 
the A in the formula is the neutron wavelength relative 
to the boundary of the slab" . 

Obviously, the neutron case differs from the case of 
atom waves only by the fact that s-wave (I = 0) scattering 
dominates the forward scattering amplitude, which, in 
this case, is given by: 

exp(i<5 )sin(5o 

/(«r) = £ 

w -a(l - ik r a) (22) 



Up to now, we have not discussed the numerical value 
of the index of refraction of gases for matter waves. A 
somewhat surprising feature is that the index of refrac- 
tion of gases for matter waves 0, [1] and the index of 
refraction of transparent gases for light has similar val- 
ues, when the gas density n t is the same. In this section, 
we calculate the value of the [n — l)/nt ratio from first 
principles in both cases and we compare the results, in 
order to understand the origin of this similarity. 

J. Schmiedmayer et al. [J] give the values of the real 
and imaginary parts of (n— 1) of several gases at T = 300 
K for sodium matter waves with a velocity v p = 1000 
m/s. All the values of lZe(n — 1) or of Xm(n — 1) are in 
the range (0.14 — 2.49) x 10~ 10 for a gas pressure equal 
to 1 mTorr at T — 300 K. This pressure corresponds to a 
gas density n t ~ 3.2 x 10 19 m~ 3 from which we find that 
the ratio lZe(n — l)/«t of Tm(n — l)/n t is in the range 
(0.4 — 8) x 10~ 30 m 3 . Our results for lithium waves Q 
are somewhat larger, in the range 1 — 2 x 10~ 29 m 3 . 

It is well known that the index of refraction of air for 
visible light is almost purely real with its value given by 
(n — 1) ~ 2.8 x 10~ 4 at ordinary pressure and a tempera- 
ture of 288 K, with a density n t ~ 2.55 x 10 25 m~ 3 . From 
these values, we calculate a ratio (n — l)/n t ~ 1.1 x 10~ 29 
m 3 for air. Similar values will be obtained for other trans- 
parent gases. 

The ratio (n — l)/n t has comparable values for light or 
matter waves, although the only common feature is that 
they involve an index of refraction. We now evaluate 
these two ratio from first principles. 
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In optics, the index of refraction of an atomic or molec- 
ular gas is dominated by the electric dipole transitions, 
with its value given by: 

n(w) = [1 + 4mi t a{uj)} 1/2 ks 1 + 2vm t a(w) (24) 

where n t is the gas density, a(tu) the atomic or molecular 
electric polarizability for an angular frequency to. In the 
visible region of the spectrum, where u> is usually smaller 
than the angular frequency of the main resonance tran- 
sitions, a(uj) is close to its static value a(0). We now use 



atomic units, with a(0) 



radius (a « 0.529 x 10 10 m) and we get: 



t (0), where ciq is the Bohr 



n(u>) — 1 



27raQa ua (0) 



(25) 



Tabulated values of a a u for atoms [32j vary from 1.4 a.u. 
for Helium up to 450 a.u. for Francium and the polariz- 
ability of small molecules has similar values (for instance, 
nitrogen dimer N 2 has a polarizability a ua (0) 11.9 

For matter waves, the index of refraction is given by 
equation (|20[) . The imaginary part Ira(n-l) is related to 
the total cross-section a(v r ) and, for a purely attractive 
— Ce/r -6 potential, there is a closed form expression of 
tr(vr) 1: 



er(iv) = 8.08 



flVr 



2/5 



(26) 



Tm{n — 1) 



= 4.12 x 10 _3 a2 



7/5 

-)(-) (C 



y 6ua ) 



,2/5 



(28) 

Tabulated values of the atom-atom Cg au coefficient span 
the range from 1.47 for Helium up to ca. 7260 a.u. for 
Cesium [33|, HH (we quote these references which pro- 
vide a large set of Cg values). The comparison of equa- 
tion (|25|) and (|28|) proves that, the index of refraction of 
gases for light and matter waves have very different ex- 
pressions and, if they have comparable values as in the 
cases discussed above, this must be considered as a pure 
coincidence. The index for matter waves has a rapid ve- 
locity dependence in v p 7 ^ 5 in a large velocity range and 
it would be considerably larger for lower projectile ve- 
locities, provided that the target gas velocity can still be 
neglected. 



IX. CONCLUSION 

In this paper, we have shown that the equations pro- 
posed by R. C. Forrey et al. [H, Ell and by our group 
[10 ] to take into account the motion of the gas particles 
in the calculation of the index of refraction for matter 
waves are not consistent with the traditional definition 
of a cross-section. Follo wing the result obtained by C. 
Champenois in her thesis [l7| , we propose a new formula 
for the index of refraction: 



The 8.08 factor is the numerical value of a complicated 
expression involving gamma function. This result is valid 
in an intermediate range of energy, with ma ny p artial 
waves contributing to the scattering amplitude [11]. This 
result does not explain the glory oscillations, which ex- 
ists only when the potential is attractive at long range 
and repulsive at short range, but this result gives the 
correct value of the cross-section averaged over the glory 
oscillations. As we want only an order of magnitude, we 
will make some simplifications, by neglecting the ther- 
mal motion of the target particle and by assuming that 
m p <C 77i t , so that we can replace v r by v p . We thus get 
an expression of the imaginary part of (n — l)/nt. 



Im(n - 1) _ 4.04 



C 6 



2/5 



1 + 27rn t - 



rni 



m t //(/., ) 



L.2 
P 



(29) 



This formula is consistent with classic results of collision 
theory and it is also in agreement with a result recently 
obtained by K. Hornberger using the formalism of the 
quantum Boltzmann equation [24j . We have used this 
formula when comparing our measurement of the index 
of refraction of gases for lithium waves with theoretical 
values 18j. Equation (|29p agrees also with the formula 
giving the index of refraction of matter for neutrons. 

Finally, in the limiting case where the thermal velocity 
a of the gas is considerably smaller than the projectile 
velocity v p , equation (1291) and equation ([5]) proposed by 
R. C. Forrey et al. [l5l . 1 1 81 ] are equivalent but they differ 



(27) considerably in the opposite limit v p <C a. 



For a purely attractive ~Cs/r~ b potential, the real and 
imaginary parts have similar values: Schmiedmayer et 
al. [J] calculated the ratio p — lZe(n — \)/Im(n — 1) and 
found that it is constant and equal to p — 0.7265 (see also 
[l6l]). We can now express Im(n — l)/n t , using atomic 
units for Cq = CQ au m e c 2 aj s ag {m e is the electron mass 
and ccfs ~ 1/137.037 is the fine-structure constant). We 
get: 
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XI. APPENDIX A 

We calculate the distribution P(v r ) of the relative ve- 
locity v r , starting from a normalized Maxwell-Boltzmann 
distribution of the target velocity v t : 



We take the z axis along to the projectile velocity v p , 
we use spherical coordinates for v r , v r being its modu- 
lus, 9 its angle with the z-axis 9 and and ip its azimuth. 
Integration over ip and 9 is easy and we get: 



PMs(v t )rf 3 v t = 



1 

7r 3 / 2 a 3 



exp (-v^) cf v t (30) 



with a — yj2ksT /mt, &b being the Boltzmann constant 
and T the temperature. Using v f = v p — vy, we can 
write: 



1 



P(v r )dv r 



2v r 



7r 1 ' 2 au 



■ exp 



«p + v r 



sinh 



2v p v r 



P(v r )d 3 v r = . exp [— (v„ — v r ) 2 /a 2 l d 3 v r (31) which is normalized, f n °° P(v r )dv r = 1. 



(32) 
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